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A WEB-portal HepWeb allows users to perform the most popular calculations in high energy 
physics - calculations of hadron-hadron, hadron-nucleus and nucleus-nucleus interaction 
cross sections as well as calculations of secondary particles characteristics in the interactions 
using Monte Carlo event generators. The List of the generators includes Dubna version of 
the intra-nuclear cascade model (CASCADE), FRITIOF model, ultra-relativistic quantum 
molecular dynamic model (UrQMD), HIJING model, and AMPT model. 

Setting up the colliding particles/nucleus properties (collision energy, mass numbers 
and charges of nuclei, impact parameters of interactions, and number of generated events) 
is realized by a WEB interface. A query is processed by a server, and results are presented 
to the user as a WEB-page. 

Short descriptions of the installed generators, the WEB interface implementation and 
the server operation are given. 



Introduction 

There are some servers for on-line calculations in particle physics and nuclear physics. The most advanced 
one is a NRV-server (Nuclear Reactions Video) [II at the Joint Institute for Nuclear Research (JINR, Dubna, 
Russia) . The NRV is a system of management and graphical representation of nuclear data and computer 
simulations of low energy nuclear dynamics. It consists of a complete and renewed nuclear database and 
well known theoretical models of low energy nuclear reactions altogether forming the "low energy nuclear 
knowledge base". The NRV solves two main problems: fast search for and visualization of experimental 
data on nuclear structure and nuclear reactions, and ondine application of commonly used models of nuclear 
dynamics. The JavaScript technology is used for forming and managing the nuclear database, and running 
some nuclear models. Thus results are directly accessible through the net with any computer having a 
Web-browser supported Java codes. 

Jet Web-server (a WWW interface and database for Monte Carlo tuning and validation) was created in 
2002 PJ. The aim of the project was "to allow rapid and reproducible comparisons to be made between 
detailed measurements at high-energy physics colliders and general physics simulation packages". The "gen- 
eral" purpose Monte-Carlo simulation programs were PYTHIA [3] and HERWIG g]. Only HERA, LEP and 
Tevatron experimental data were included. 

At the centre of JetWeb there was a JAVA object model containing the properties and interactions of 
Models, Papers, Plots and Fits. The data underlying this model were stored in the JetWeb database - a 
MySQL database. A "Model" completely specified a unique generator, version and set of parameters. A 
"Paper" encapsulated the measured data from a single publication and was associated with measured cross 
sections. A "Fit" contained the results of a comparison between real data and predictions of a specific model. 
The data and MC results were processed with the help of HZTOOL package [5. , and output histograms in 
the form of XML files in the Java Analysis Studio (JAS) [6] plotML DTD format [7| were written. The 
plotML data were converted into the Java object model, and written to the JetWeb database via JDBC [S]. 
The user accessed JetWeb via the web interface, which consisted of Java servlets run on a Tomcat [5] server, 
delivering HTML pages written using the JetWeb HTMLWriter facility (ucl.hep.jetweb.html package). The 
servlets accessed the JetWeb database via the JDBC calls encapsulated within the java object model. 

If a fit requested by the user were done before, the static webpages were searched for and send to the 
user. If there were no results stored for the user's specified set of parameters, a new job request with the 
required parameters was generated. 

Now JetWeb does not operate. It seems to us that the matter is that the HZTOOL library of the 
experimental data and corresponding FORTRAN routines has not been update for a long time. As the data 
set was restricted, all calculations were performed quite fast, and there was no need and possibilities to make 
new ones. 
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We think that the main idea of the Jet Web server was quite interesting, and we use them at the creation 
of our server the Hep Web server. 



The other interesting server is "Qt resummation portal at Michigan State University (http://hep.pa 
msu.edu/wwwlegacy/) created at 2003 [TO]. It allows one to plot transverse momentum distributions for 
cross sections of several particle reactions. The following processes are implemented there: 

• Massive vector boson production - pp — > W X, pp — > Z°X; 

• Photon pair production - pp — > 77X; 

• Z-boson pair production - pp — > Z°Z°X; 

• SM Higgs boson production - pp — > H°X. 

The output figure shows distributions da / dQ 2 dydqx for the production of on-shell particles (or pairs of 
on-shell particles in the case of the 77 and ZZ production) with specified invariant mass Q, rapidity y and 
transverse momentum qx in the lab frame (the center-of-mass frame of the hadron beams). A user can 
obtain resummcd, fixed-order and asymptotic cross sections. Details of program implementation is unknown 
for us. 

In Sec. 1 below we consider tasks of on-line calculations and formulate the aim of our project. A 
description of how to use our server for Monte Carlo calculations in CASCADE, FRITIOF, UrQMD and 
HIJING models is given in Sec. 2. Technical details of the portal implementation are considered in Sec. 3. 
Appendices present short descriptions of the installed generators. 



1 Task of on-line calculations 

Monte Carlo event generators play a very important role in high energy physics. One can mark the following 
areas of their application: 

1. Pragmatic or practical tasks - development of new or upgrade of old experimental setups to study 
some processes/interactions, design of detectors, Monte Carlo simulation of the detector responses and 
so on. Event generators applied at these should have fast operation speed, stability of work, and a 
rough reproduction of previous experimental results. As an example, let us mention applications of the 
UrQMD model [TT1 IT2] for the development of detectors for research on nucleus-nucleus interactions 
(CBM collaboration [13]), and detectors for investigation of antiproton-proton annihilations (PANDA 
collaboration [H]) at future GSI accelerators. The RQMD [TS] and HIJING QUE] models have been 
used for analogous purposes for RHIC experiments. The well known Geant4 package |18j is widely 
used to simulate various installations. 

2. Analysis of new experimental data and planning of new investigations . They include 
comparison of new data with previous data and model predictions. As a rule, the new data do not 
agree with model predictions, so some questions arise in this regard: whether all special features of the 
setup have been taken into account; whether they are free of methodical errors; whether the theoretical 
models are used correctly; whether the model parameters set is right; whether the discrepancy between 
the experimental data and model predictions is of systematical character; whether the discrepancy was 
observed in previous experiments; whether the discrepancy was considered as an evidence of a new 
physical effect, and so on. They are solved differently and often quite difficult using quite often event 
generators. The generators should have a flexibility in parameters variation and physical scenario, as 
well as a sufficient physical meaning of the parameters. 

Another situation takes place at new research planning. A first question asked by experimentalists 
is related to the load of the setup by ordinary background processes. A second deals with the radiation 
condition of the experiment. A third question asks about the admixture of the background processes in 
the phenomenon under study, and how it can be damped, and so on. Clearly, experimentalists prefer 
to use well approbated and well recommended models for their estimations. Here one can not avoid a 
study of model application experience. As a rule, there is no enough time to do it. Thus, experimental 
collaborations attract the authors of the models or use authors variants of model code to solve the 
questions. As an example, let us point out on the estimations of secondary particles multiplicity in 
central gold-gold interactions obtained by CBM collaboration [19] , and the estimations of background 
processes intensities by the PANDA collaboration [2"U] . 
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3. Scientific or cognition aims only - search for new effects or phenomena on the base of analysis of 
a discrepancy between experimental data and model predictions. One uses the fact that Monte Carlo 
models are a synthesis of existing notions about process mechanics. Thus, the discovered discrepancy 
can be considered as an evidence of our insufficient understanding or else as an evidence of new effects. 
For example, the discrepancy between experimental data and intra-nuclear cascade model calculations 
growing with the collision energy rising led in its time to appearing a very important conception for 
high energy physics - "formation time of secondary particle". 

The final aim of all the efforts is creation of a theory of processes that could predict effects with any 
predetermined exactness. As there are only few such theories, the aim is re-formulated - creation of a theory 
or a model predicting observable effects with specified exactness. Determination of the exactness is a special 
additional task. 

The philosophical aspect of the scientific research - the " cognition of Good wisdom" , is out of the scope 
of our consideration. 

The aim of this work is to create a WEB-portal which allows one to perform calculations in the second 
stream. 



■T^pr 

hepweb" 




Monte Carlo Simulation Tools for High Energy Physics 

The page allows you to perform Monte Carlo simulations of high energy physics processes, and to evaluate the main properties of interactions: hadron-hadron, had r on -nucleus and nucleus -nucleus interactions. 
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Figure 1: View of the main page. 



Though the variety of tasks of the second stream is quite bright one can select out some common tasks. We 
think, a user can be interested first of all in global properties of interactions such as interaction cross sections, 
charged particles multiplicity distributions, momentum, angular and energy distributions of particles. There 
is a quite defined need in theoretical calculations of impact parameter distributions, multiplicity distributions 
of intra-nuclear collisions, and multiplicities of wounded nucleons. It is desirable for check of a model quality 
to know whether energy, momentum, baryon number, lepton number and so on are saved in the model or 
not. 

At high energies total, inelastic and diffraction cross sections are usually calculated in the Glauber 
approach. The reggeon theory can be used also for the same aim. All of the possibilities are offered by 
our server. There is a possibility in the Glauber approximation to calculate impact parameter distributions, 
multiplicity distributions of intra-nuclear collisions, and multiplicities of wounded nucleons. 

The most easiest way to calculate the inclusive and global properties of interactions is utilization of event 
generator programs. There are a lot of event generators especially at low and intermediate energies which 
are of great practical interests. Most of them arc not accessible. It is impossible to review the left part of 
them. At high energies a set of the generators is not so large. A list of most actual ones is quite restricted. 
At energies below 10 GeV/nucleon various variants of intra-nuclear cascade model are very popular. We 
have installed at out server the variant of CASCADE code [21] quite well known in JINR. At high energies 
the FRITIOF model [23] was applied quite well in the past. At present time the ultra-relativistic quantum 
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molecular dynamics model (UrQMD [TT1I12) ) is applied very often. At super high energies, the HIJING model 
QUE], PYTHIA [3] and HERWIG 0] models are utilized. The PYTHIA and HERWIG can be applied only 
for simulation of elementary particles interactions. In order to simulate nucleus-nucleus interactions most 
actual now we have installed at our server the FRITIOF, UrQMD and HIJING models. The other models 
will be included in future whether according to user wishes, or mastering of new codes. 

The models selected by us are quite well documented for experienced user as a rule. Though we think 
that most of the potential users use the model episodically. For them we have created a simplified access to 
the generators using WEB-based interface. 

The main page of the portal looks like it is presented in Fig. 1. The address in INTERNET is: 
|http://hepweb.jinr.ru/| 

2 User guide 

2.1 Cascade Evaporation Model(CASCADE) 

The page presented in Fig. 2 gives a user a possibility to calculate properties of interactions using well known 
cascade-evaporation model [23 - 28J. A short description of the model is provided by the page also. The user 
should set up mass numbers and charges of colliding nuclei, an interval for the impact parameter sampling, 
a momentum of a projectile nucleus per nucleon in the target nucleus rest frame, a number of generated 
events, and run the calculation process pushing "Get Results" button. All of these are prompted by an 
interface presented in Fig. 2. 

CASCADE - intra nuclear cascade-evaporation model 

The page allows you to perform Monte Carlo simulations of hadron-hadron, hadron-nucleus and nucleus-nucleus interactions. 



Choose the projectile and target 





Projectile 


Nucleus T 


Mass number 4 


Charge 2 


Target 


Nucleus T 


Mass number 12 


Charge 3 


Arb trsr/ rp=ct p=r=reter T 


Bmin 0. 


Bmas 10. 


Reference frame Lab. system (target at rest) 


PlabinGeV 4.2 



The code tan rtgj amsi i:i- bir^ti-hadtpn interagionsl 



Number of generated events too 

[ Get HESiiits~~) 
Short description of the model fpdf) 
If vol have any comments or suggestions please contact us at uzhi@cv.Jinr.ni 

Figure 2: View of the CASCADE page. 

If the user is going to study hadron-nucleus interactions, he/she has to choose in the "Projectile" list 
a needed item. The user should obligatory set up a mass number and a charge of the hadron. The mass 
number is treated in the case as a number of particles in a projectile. Thus it must be equal to one. 

If the user is going to study minimal bias interactions, he/she can use the following consideration for a 
choosing of the impact parameter interval. The radius of a nucleus with mass number A is near to the value 
Ra ~ 1.2 A 1 / 3 (fm). Thus for interactions of nuclei with mass numbers P and T the maximum impact 
parameter can be estimated as Rp + Rt + 2 (fm). "2" represents a double radius of NN- interactions. The 
minimal value of the impact parameter is obviously equal in the case to 0. 

For a projectile hadron the user can set up Rp to 1 (fm) (nucleon radius). 

If you know experimental cross section of the interactions, you can estimate the maximum value of the 
impact parameter as U a exp /Tr (fm). 
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In other cases you should undertake additional consideration. A Glauber cross section calculations can 
be useful for you at this (see below). 

If you set up a small number of events, and are interesting by light nuclei interactions, the results can be 
obtained quite fast. They will be presented on the page described below. 

Results:: 
Host: 

Root File 

E-nerrv ;c ==-l~ - zisi Ion 
Fx momeiicuin conservation 
Pv rrh rym gti ti_m* conserv-ation 
Vr~ momentum conservation 
Charge conservation 
Barvon number -conservation 

■A^irnuttial isotropv l"az:imutnal angle distribution) 

Impact pafimscgf distribution 

Ctiarged particle multiplicity distribution 

Charged particle pseudo-rapidin r distribution 

Charged particle rapidity distribution 

Charged particle ]Pt distribution 

Ctiarged particle energy' distribution 

Ctiarged particle CosiXheta;) distribution 

Cliarged particle IPni distribution 

IParticle compisition [IDs^ 

Evap. 5?- IN and Frag. Idn. £ distribution 

E-vap. _P. IX and Frag. CosfT^ distribution 

Figure 3: View of the Results page. 

If your calculations require a lot of time, the view of the page does not change immediately. This means 
that your task is put in a queue for execution, or is executing now. In the following you have two possibilities: 
to be at the page refreshing it periodically until the results will be ready, or left the page and turn to your 
busy. 

The results of your task execution will be stored in a data-base. To reach an access to them in a new 
connection you should set up one more just the same parameters of the task (mass numbers and charges of 
nuclei, impact parameter interval, momentum of the projectile, number of events), and should run process 
as it was described before. The server looks first of all though the data-base. If the corresponding results 
are in the data-base, they will be presented on a new pag^J . 

The calculation results are presented on a page, a view of which is given by Fig. 3. In the upper part 
of the page, interaction parameters given by a user are presented. After that a list of calculated interaction 
characteristics are followed. At clicking on an item of the list an additional window will be opened with a 
graphical representation of a characteristic. The page give you a possibility to download the results as a 
ROOT files. 

The version of the cascade-evaporation model developed by Zh.Zh. Musulmanbekov [5T] is used at the 
calculations. A short description of the model extracted from Ref. can be found at the CASCADE page. 

2.2 FRITIOF Model 

A page presented in Fig. 4 gives you a possibility to calculate properties of interactions using the famous 
FRITIOF model (25]. The code of the model was essential modified by V.V.Uzhinsky [301: calculations of 
cross sections and a sampling of interacting nucleons within the Glauber approximation were added; an unit 
for a simulation of the nuclear destruction at the fast stage of interactions (the reggeon nuclear destruction 
model[31J) was added also; the evaporation of nuclear residues was taken into account and so on. Most of 
the changes are described in Ref. [2J5]. An extraction from the paper related to the model is presented 
on the page. 

To run the simulation process you should set up: 

1 ' Unlike the other portals we do not undertake any efforts to identify users and their tasks. We are thinking that all 
calculations are unique ones, and their results are important. 
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Modified FRITIOF 

The page iIIqtfs you to perform Monte Carlo simulations of hadron-hadron, hadron-nucleus and nucleus-nucleus interactions. 

Choose Che projectile and target 



[Projectile | Nudeus ▼ | Mass number 4 




| Charge 2 




[Target | Nucleus T | Mass number 12 




| Charge G 






| Arbitrary impact parameter T | EfQtQ 0. 


| BmsK 10. 








Nuclear destruction parameter- strong (1) or week (0) 






| CND 0-2 


Fsrffli momentum in GeV c 






0.2 



ONLY STABLE PARTICLES RECORDED 



[Reference frame Lab. system ▼ [piab or Sqrt<s) in GeV 4.2 



Number of generated events 100 

[ Get Results ] 
S kort des cription of the ni o del ( pdfl 
If you have any comments or suggestions please contact us at uzhi@cv.jirr.ru 

Figure 4: View of the FRITIOF page. 



• Mass numbers and charges of colliding nuclei, or choose a projectile particle from a pop-up menu (in 
the case arbitrary mass and charge are needed also); 

• If you are interested in the minimal bias interactions miss the next line. In other case, choose from a 
pop- up menu the item - "Impact parameter in the range", and set up minimal and maximal values of 
the interaction impact parameter. They must be meaningful, see the previous section. 

• We do not recommend to change the item - " Nuclear destruction parameter" , without a special need. 
It allows to variate the nuclear destruction power at the fast stage of interactions. The prompted value 
is determined at an analysis of heavy nuclei interactions. For light nuclei interactions it can be changed 
to 1. 

• The item - "Fermi momentum", allows you to fit the Fermi- momentum for a description, for example, 
the spectator fragment spectra. 

• The following pop-up menu gives you a possibility to order what to do with unstable particles at the 
end of the simulation of the fast stage of interactions. 

• Particle characteristics can be calculated in the laboratory system, or in the centre-of-mass system. De- 
pending on your choice, set up a momentum per nucleon of the projectile nucleus, or the corresponding 
energy of the TV TV-interactions. 

• It is natural that you should set up a number of the generated events. 

• At the end of your typing, click on the button - " Get Results" . 
Calculated characteristics will be presented as it was described before. 



2.3 The Ultra-relativistic Quantum Molecular Dynamics Model (UrQMD) 

The page presented in Fig. 5, gives you an opportunity to calculate the interaction characteristics using 
UrQMD model [TJJ [TJ]. Changes, made in the original code of the model, are described in details in 
Ref. [32]. 

To run the simulation process you should set up: 

• Mass numbers and charges of colliding nuclei, or choose a projectile particle from a pop-up menu (in 
the case arbitrary mass and charge are needed also); 
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IheUtrarelalivislic Quantum Molecular Dynamics model is a microscopic model used to simulate (ultra)relatiristic heavy ion collisons in the energy range from Bevalac and SIS up to AGS. SPS andRHIC, 




The page allows you to perform Monte Carlo simulations of hadron-hadron, hadron-nucleus and nucleus-nucleus interactions. 
Choose the projectile and target 



Projectile ^ T 


Mass number 1 


Chaise 1 


Target Hxkus w 


Mas* number 1 


Charge 1 


Arbitrary impad parameter * 


Rmin 0. 


Bmax io. 




|Time of evalutiaa tn fm'c IW 










|Time of output in fm'c ,M 








(Type of tateractions.: 


Inelastic collisions without potential interactions t 









Point the collision energy 



ference frame Lab. system " (plaborSqrE(s)iii<jeV 4-2 

Nomberofgeneiatedeventj 1m 



Figure 5: View of the UrQMD page. 



• If you are interested in the minimal bias interactions miss the next line. In other case, choose from a 
pop- up menu the item - "Impact parameter in the range", and set up minimal and maximal values of 
the interaction impact parameter. They must be meaningful, see the previous section. 

• We do not recommend to change the items - "Time of evolution in fm/c", and "Time of output in 
fm/c" without defined need. They determine a time of interactions and a frequency of an intermediate 
results output. It is obvious that the time of evolution must not be lower than the time of a particle 
penetration through a target nucleus. At low and intermediate energies the change of the parameter 
has an influence on properties of the pre-equilibrium particles and residual nuclei. Giving a small time 
of the evolution, you will obtain a "hot" residual nucleus with a large mass because the pre-equilibrium 
particles will be captured by the residual. At large evolution time the nucleus will be more "cold", and 
the multiplicity of the pre-equilibrium particles will be larger. 

• The item - " Type of interactions" , allows you to order to perform calculations with or without potential 
interactions. The corresponding sub- menu items are "Inelastic collisions + potential interactions" 
and "Inelastic collisions without potential interactions". At energies larger then 4 GeV/nucleon all 
calculations are performed without the potential interactions. 

• Particle characteristics can be calculated in the laboratory system, or in the centre-of-mass system. De- 
pending on your choice, set up a momentum per nucleon of the projectile nucleus, or the corresponding 
energy of the iViV-interactions. 

• It is natural that you should set up a number of the generated events. 

• At the end of your typing, click on the button - " Get Results" . 
Calculated characteristics will be presented as it was described before. 

Below the menu there is a reference on the UrQMD model validation WEB-page which collects materials 
on a comparison of the UrQMD model calculations and experimental data [33 . 



7 



2.4 HIJING Model 



The page presented in Fig. 6, gives you an opportunity to calculate properties of hadron-hadron, hadron- 
nucleus, and nucleus-nucleus interactions at high and super high energies using the HIJING model [17]. We 
use a variant of the code disposed at the code depository of the LCG project [34] and described in Ref. [35] . 



HIJING - A Monte Carlo Program for parton and particle production in high energy hadronic and nuclear collisions 

The page allows you to perform Monte Carlo simulations of hadron-hadron, hidron nucleus and nucleus -nucleus interactions, 

Choose the projectile and target 



Projectile Nuc^ie - 


Mass number i 




| Charge 1 


Target Nucfee 


Mass. number i 
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Impact parameter Min o. 
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[ 



P.e:eren:r xams CMS ■> 



Change parameters of the model, if it is needed, 
List of the parameters is given below. 
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HI JL\G validation ho me -pag e 



Figure 6: View of the HIJING page. 



To run the simulation process you should set up: 

• Mass numbers and charges of colliding nuclei, or choose a projectile particle from a pop-up menu (in 
the case arbitrary mass and charge are needed also); 

• If you are interested in the minimal bias interactions miss the next line. In other case, choose from a 
pop- up menu the item - "Impact parameter in the range", and set up minimal and maximal values of 
the interaction impact parameter. They must be meaningful, see the previous section. 



Set up a number of generated events. 



• We do not recommend to change one of the most important parameters of the model 
transverse momentum of hard scatters", if you have not a special need. 



' Minimum 



• Particle characteristics can be calculated in the laboratory system, or in the centre-of-mass system. De- 
pending on your choice, set up a momentum per nucleon of the projectile nucleus, or the corresponding 
energy of the TV TV-interactions. 

• We give in the menu an additional possibility to change the model parameters and options. A list of 
the parameters and options is given below. For example, to switch off the jet quenching effect it is 
need to put IHPR2(5)=0. This means that you should enter the corresponding windows numbers "5" 
and"0". 

• At the end of your typing, click on the button - " Get Results" . 
Calculated characteristics will be presented as it was described earlier. 

Below the button - " Get Results" , there is a reference on the HIJING model validation WEB-page [55] 
which collects materials on a comparison of the HIJING model calculations and experimental data. 
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2.5 Glauber cross section calculation 

A page presented in Fig. 7 gives you an opportunity to calculate hadron-nucleus and nucleus-nucleus inter- 
action characteristics in the Glauber approximation such as the total cross section, the elastic cross section, 
the AB— >XB cross section, the AB— >AX cross section, the AB— >X cross section, the production cross 
section, the production cross section of interactions with composite materials, and various distributions of 
inelastic reactions. For this you should choose needed item and click a button - "Run". 

Glauber's Cross Sections 

The page allows you to evaluate Glauber 1 s eras s s ections for hadron-nucleus and 

nucleus-nucleus interactions 







Total cross-section 


Run! 




Elastic cross-section 


Run! 


AB->XB cross-section 


R-t! 


AB->AX cross-ssction 


r_-: 


AB->X cross-ssction 


Run! 


Production cross-ssction 


Run! 


All above listed cross sections 


Run! 


Production cross-section of int tactions with composite materials 


Run! 


Various distributions of inelastic reactions 


Run! 



Choose needed, and push the corresponding button 



Figure 7: View of the Glauber page. 

All following pages for the Glauber calculations have nearly the same view presented in Fig. 8. To perform 
a calculations you should set up: 

• Mass numbers and charges of colliding nuclei. If you are going to study hadron-nucleus interactions, 
set up for the mass and charge of the projectile the values " 1" and " 1" . 

• After that you should set up AW-interaction properties - the total interaction cross section, the slope 
of differential elastic scattering cross section, and the ratio of real and imaginary parts of the elastic 
scattering amplitude at zero momentum transfer. All of the quantities can be found in the compilation 
|36j . The same quantities can be calculated assuming the gaussian parameterization of the elastic 
scattering amplitude and ReF/ImF — if you know the total and elastic NN interaction cross 
sections. For a determination of the NN cross sections one can use parameterizations proposed by the 
Particle Data Group (PDG [37] )■ We allow both of the possibilities. To choose one of the methods 
you should mark one of them in the corresponding check box. 

If you choose an usage of the PDG paramerization, you should additional set up a reference frame and 
a momentum per nucleon or an energy in the centre-of-mass system (y/s). 

You need to set up the required parameters by yourself if you are going to consider hadron-nucleus 
interactions (it A, KA). 
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• In the codes the stochastic method of an evaluation of the multi-dimensional integrals is used. It is 
quite well described in Refs. (101 FfTl W2\ . Thus, you have to set a number of samples (statistic). It 
determines an accuracy of the calculation. 

• At the end of your typing, click on the button - " Get Results" . 

Glauber's Cross Sections 

All cross sectiaas 

[Enter projectile mass number and charge: |a=1 [^1 

(Enter target mass number ffldcharje: |e=i |?i 

You can use Particle Data Group paramet erization 
of NN cross-sections, or introduce them by yourself. 

Choose one of them. 



PD G parameterization 


□ 


You can use CMS energy or Plab: 


Eons 9 100. GeV 



|My parameters 


□ 


|NN total cross-section: 


40. (mb) 


[Slope of elastic cross section at t=0 
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Figure 8: View of the Glauber calculation page. 



2.6 Reggeon cross section calculation 

The page presented in Fig. 9 gives you an possibility to calculate hadron-hadron, hadron- nucleus, and 
nucleus-nucleus interaction properties using the reggeon or pomeron approach. 
To perform a calculations you should set up: 

• Mass numbers and charges of colliding nuclei. If you are going to study hadron-nucleus interactions, 
set up for the mass and charge of the projectile the values "1" and "1". If you are interesting in 
nucleon-nucleon collisions set up also for the target the values "1" and "1" (Proton is a nucleus of 
hydrogen atom.) 

• You should set up a reference frame and a momentum per nucleon, or an energy of iVW-interactions 
in the centre-of-mass system. They are needed for a calculation of the reggeon parameters. 

• Because in the code the stochastic method of an evaluation of the multi-dimensional integrals is used, 
you have to set a number of samples (statistic). It determines an accuracy of the calculation. 

• At the end of your typing, click on the button - " Get Results" . 

A description of the calculating method will be published elsewhere. 
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Reggeon Cross Sections 



Enter projectile mass number and charge: 


A=4 


ZP2 


Enter target mass number and charge: 


E=4 


Z=2 



Point the collision energy 



Reference frame CMS ■ 



Collision energy !40W (GeV) 



Point the statistic 



Inter statistic (500 - 1000) 



In the cods stochastic method of evaluation of multi-dimensional integral is used. 
Thus, you have to set number of samples (statistic). It determines accuracy 7 of the 
calculation. 



Figure 9: View of the Reggeon page. 



3 Structure of WEB-service Hep Web 

Already existing servers described in the Introduction have been analyzed, and the following requirements 
to the Hep Web service have been formulated: 

1. The Hep Web service should have a WEB-oriented architecture; 

2. Requirements to a user's computer should be minimal one both on hardware and software levels; 

3. Generators should have standardized input and output streams; 

4. Generators satisfying to the point 3 should be easily installed and modified on the server, i.e. they 
should not have any special settings necessary only for the Hep Web; 

5. Result of a generator work should be presented in standard graphic formats with possibility to create 
a root file or text fragments; 

6. Result of the generator work should be saved because some calculations require a lot of a computational 
time, and it is inexpedient to repeat them. 

The following structure was accepted and implemented in the accordance with the given above require- 
ments: 

• The Hep Web service is based on the "client - server" architecture. 

• The HepWEb service is implemented in Java, and uses the JavaServer Pages (JSP) technology. This 
allows one to result as XHTML - or XML-documents. User can apply any Web-browser (IE, Mozilla, 
Opera etc.) for their viewing. 

• The output streams of all Hep Web generators are standardized. 

An output data is written in a text file. Any number of histograms can be stored in the file. A 
histogram has a header where a name of the histogram and the names of the co-ordinate axes are 
given. Then a set of co-ordinates is followed. Additionally a ROOT-file of the histograms could be 
also created. 

• The results are stored in the database. 
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All generators are presented on the server as executable modules (the source files of the codes are 
unavailable for users). The modules can be also started independently of the HEP WEB. This allows one 
to upgrade generators without reboot of the server and add new generators specifying only ways to the 
corresponding executable modules in a configuration file. 

The standardized output stream is analyzed by a special method created by us in order to extract the 
histograms. The histograms are converted into graphic files (png-files) using the standard Java-technology. 

The MySQL database was chosen to store the results because it is free, and its bandwidth of the read/write 
data is enough for the HepWeb service. 

A interaction scheme between the components of the HepWeb service in presented in Fig. 10. 
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Figure 10: Structure of the Hep Web-server. 



A client generates a request to the server using a form which is built in a static HTML page. Each 
generator uses its own HTML page with an according form. An answer on the query is generated by the 
server with a help of the JSP technology. The server checks correctness of the parameters of the query. For 
example, it is checked that the energy of iV./V-collision in the centre-of-mass system y/s is larger than the 
sum of nucleon masses, or that the range of the impact parameter does not extend to negative values, etc. 
All correct requests are stored in the database. The server checks for each request if a request is already 
stored in database or not. If the request is stored, then the server generates the result page and sends it to 
the user. If the results for the request are not presented, an empty page is sent to the client with a line - 
" Please wait a few minutes then refresh this page" . In the case the requested generator with corresponding 
parameters is started on the server. 

The HepWeb service has two possibilities for executing the generator: 

1. Use only one computer for the server; 

2. Use distributed computing system for the server. 

3.1 The single-computer based HepWeb service structure 

For the structure, all calculations are performed by one computer. A request to run a generator follows 
to the consecutive list. Thus, all queries of a generator will be executed in consecutive order. The given 
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structure is efficient for a small number of the requests. If the number of queries increases, the astronomical 
time of their execution will increase too. The structure is marked by letter "A" on Fig. 10. 

3.2 Structure of the Hep Web service based on a distributed computer system 

The structure is developed for starting tasks on several computers. It allows one to perform several queries 
simultaneously. In the presence of a sufficient quantity of computers in the system, the performance time of 
a new query is nearly equal to the performance time of this task on one computer. This structure is marked 
by letter "B" on Fig. 10. 

The Hep Web service was tested in 2006 on distributed metaclaster of the Dubna-Grid project ( [38] . [39]). 
An interaction of the service in batch-mode of access to resources of the metaclaster has been successfully 
proved. 

Summary 

A WEB-portal is designed which allows a user to perform the most popular calculations in high energy 
physics - calculations of hadron-hadron, hadron-nucleus and nucleus- nucleus interaction cross sections and 
calculations of secondary particles characteristics in the interactions using Monte Carlo event generators. 
The list of the generators includes the intra-nuclear cascade model (CASCADE), the FRITIOF model, the 
ultra-relativistic quantum molecular dynamic model (UrQMD), and HIJING model. 

Setting up of colliding particles properties (collision energy, mass numbers and charges of nuclei, impact 
parameters of interactions, and a number of generated events) is realized by a WEB interface. A query is 
processed by server, and results are presented to the user as a WEB-page. 

Appendix: The Cascade-Evaporation Model (CEM) 

The cascade-evaporation model (see Refs. [23]-[2S]) was one of the first models of hadron-nucleus interactions 
at high energies (E > 1 GeV) very important for practical applications. It saves its position at present time 
too. The model assumes that due to an inelastic interaction of a projectile hadron with one of target nucleons 
a new particle is produced. The participating target nucleon accepts momentum and begins moving in the 
nucleus. All moving (cascade) particles can interact with other nuclear nucleons to produce new particles 
or suffer an elastic re-scattering. Therefore, a cascade reproduction of moving particles is assumed. The 
interactions between the cascade particles are usually omitted. The process continues until all moving 
particles either leave the nucleus or are absorbed. 

In the case of nucleus-nucleus collisions, it is assumed that the cascade particles can interact with 
projectile and target nucleons. To choose the nucleons which can participate in elementary interactions, 
all of the nucleons of the colliding nuclei with mass numbers A and B are identified by the coordinates 
{xi,yi, Zi, 1 < i < A) and (x'p y'-, z'j, 1 < j < B) in the initial state in the corresponding reference frames. 
At the sampling of the nucleon coordinates for a nucleus with A (B) < 10 the oscillator density is used. 
The Wood-Saxon density with the parameters Ra = 1-07 * A 1 / 3 fm and c = 0.545 fm is applied for more 
heavy nuclei. The nucleon core is taken into account at sampling the nucleon coordinates according the the 
densities pA and ps (two nucleons can not be be closer than R c , R c = 0.4 fm). 

Taking into account the Lorentz contraction of the projectile (A) in the rest frame of the target nucleus, 
the corresponding coordinates are redefined as Zi —^Zi/'j — Ra/i ~ Rb- Here 7 is the Lorentz factor of the 
projectile nucleus, Ra and Rb) are radii of the nuclei. If the coordinates of the i-th nucleon of the nucleus 
A and the j-th nucleon of the nucleus B satisfy the condition 

(b x + Xl - x' ) 2 + (b y + yi - y'j) 2 < (R mt + A D ) 2 , 

then the two nucleons are considered as possible participants in the elementary interactions. Here (b x , b y ) 
are components of the impact parameter, Ri n t is the strong interaction radius (1.3 fm) and Xd is the de 
Broglie wavelength of the projectile nucleon. The possible participant can fly without interaction, or suffer 
an elastic scattering, or an inelastic interaction. The corresponding probabilities are 

W tot (ir * (R int + \ D ) 2 - a tot )/(ir * (R mt + X D ) 2 ), 

W el =a el /(Tr*(R mt + \ D ) 2 ) 1 
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W in = a in /(7r*(R int + X D f), 

where a tot , <r el and a m are the total, elastic and inelastic nucleon-nucleon (NN) interaction cross-sections. 
The time of possible interactions in the frame of the target nucleus is determined as t,j = (z'j — Zi)/v where 
v is the velocity of the projectile nucleus. 

A minimal time interaction will be caused first. The nucleons taking part in the interaction are ascribed 
a Fermi momentum in the corresponding reference frame and the momentum of the projectile nucleon is 
transferred to the target nucleus rest frame. It is assumed that the Fermi momentum distribution has the 
form 

W(p)dp = Zp 2 / P F {r)dp, 0<p<P F (r), (1) 
where Pf{t) depends on the local density of the nucleus p{r), 

P F = ft(37r 2 p(r)) 1/3 . (2) 

r is the radius-vector of the interacting nucleon in the corresponding nucleus. In the following the Fermi 
momenta are ascribed only to new involved nucleons. 

An elementary interaction (elastic scattering or inelastic collisions) will be rejected, if at least one of the 
two colliding nucleons, after the interaction, falls into the region of occupied states of nucleus A or B (has 
an energy lower than the Fermi energy of P F /2m). In this case the next possible interaction in the time is 
considered. This repeats until an interaction occurs. 

The simulation procedure of elastic and inelastic nucleon-nucleon and meson-nucleon interactions is 
described in detail in [23] • It allows one to reproduce the experimental observations up to energies ~ 10 - 
15 GeV, and especially the particle momentum distributions. 

After the first real interaction, the time is increased by ty. The coordinates of the moving particles 
are changed as Zk — > Zk + vtij . The coordinates of newly produced particles (tt- mesons) in the interaction 
with equal probabilities are identified by the coordinates of the target or projectile nucleons (r^ or r^). 
Then all possible interactions of the produced and moving particles with the projectile or target nucleons 
are considered. Each cascade particle can interact with nucleons in a tube with radius Ri n t + Ad along 
its trajectory. The time for possible interactions is determined as = (vk(ji — fl))/ \ Vk I or tjk ~ 
(vk(r l j — f"k)/ I I) where Vk is the velocity of the cascade particle k, and r'l is its radius vector. 

Among all possible interactions the one with the minimum time is chosen. The processes continue until 
all possible interactions have been considered. 

As mentioned above, the nucleon taking part in an interaction is considered as a cascade particle. Due 
to of this, the local density of the nucleus gets lower. This is the so-called trailing effect. 

At the end of the fast cascade stage of the process, the number and charges of spectator nucleons as well 
as the charges of the absorbed mesons determine the nuclear residual mass number and charge. 

The excitation energy of the nuclear residual is the sum of the energies of absorbed particles and the 
holes counted from the Fermi energy. It is assumed that due to an interaction of a cascade particle or a 
projectile nucleon with a target nucleon having the energy p 2 /2mjq , where p is determined according to the 
distribution |l]), a hole with an energy E^ = P F /2mN — p 2 /2m^ is created in the target nucleus. ( P F is 
given by equation ([2])). If an energy of the target nucleon after the interaction - T, is larger than the Fermi 
energy (P F /2m^) but lower than P F /2m^ + Eb (Ej, = 7 MeV), it is assumed that the target nucleon will 
not leave the nucleus. It is considered to be an absorbed nucleon. In this case its yield to the excitation 
energy of the nucleus is determined as E n = T — P F /2mpf. A meson having an energy lower than 25 MeV 
counted from the bottom of the potential well is considered as to be an absorbed one, too. Its yield to the 
excitation energy of the nucleus is determined as + m v . (It is assumed that the elementary interaction 
is caused in the potential well with a depth P f /2ttin + -E*. This value is added to the kinetic energy of 
the incoming cascade particle or the projectile nucleon and is subtracted from the kinetic energy of each 
outgoing particle.) 

The total excitation energy of the target nucleus is a sum of the energies of the holes, the energies of 
absorbed nucleons and mesons. The analogous procedure can be applied to determine the excitation energy 
of the projectile nucleus. It is obvious that the excitation energy of the target nucleus will be proportional 
to the number of ejected nucleons if one neglects the absorption of the nucleons and mesons. 

Let us note that according to the above prescription in the limit case when all nucleons of the target 
nucleus are ejected one can obtain a nucleus without any nucleons, but with a defined excitation energy. 

The excitation energy governs the nuclear residual relaxation. So the method of the excitation energy 
calculation links the fast and slow stages of the interaction. 
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The nuclear residual relaxes before thermodynamic equilibrium and can emit the so-called pre-equilibrium 
particles. This process takes place if the number of quasi-particles - N q = Nh + N n , is larger than the 
equilibrium value N q ^ eq ^ — 2y // 6aE* /tt 2 . Here Nh is the number of holes (the number of participating 
nucleons), N n is the number of absorbed cascade nucleons, E* is the excitation energy of the nucleus, and 
a is the level density parameter, taken as = A/10 MeV -1 . The pre-equilibrium decay of the nucleus is 
simulated in the framework of the cxciton model [43 . The decay of a thcrmalized nucleus is described by 
the usual evaporation approach [331 [35] . 

Additional details of the model can be found in [55]. We have installed at the Hep Web server a code 
described in [ST] , 

From the view point of modern approaches, the realization of CEM simulation is regarded as being too 
simple. It fails to take into account many important effects: a variation of the average nuclear field during the 
collision, a production of meson and baryon resonances, a finite formation time of a particle, a coalescence 
of nucleons, a multi-fragmentation of nuclei and so on. Thought the analogous implementation of the model 
[46] was recognized as the best model code used in the physics of the intermediate energies. 

Appendix: The UrQMD Model 

A variation of the average nuclear mean field and its influence on the inclusive distributions at low and 
intermediate energies are commonly considered in the VUU/BUU approaches (Vlasov-Weling-Ulenbek or 
Bolzman-Weling-Ulenbek, see Refs. [48j|49]). At higher energies or at a strong breakup of nuclei the effects 
of the average field render weak. Here, the quantum molecular dynamical approach considering the explicit 
form of the two- and three particle interactions is preferably applied. Its relativistic generalization combined 
with quark ideas of the multi-particle production is presented in the RQMD model [5TJ [S3] and UrQMD 
model [HHH]). It allows one to describe numerous characteristics of produced particles in nucleus- nucleus 
(AA) collisions at high energies. The UrQMD model is well described in [TTJ [12]. Thus below we give only 
the main features of the model. 

The UrQMD model describes hadronic interactions at low and intermediate energies (-y/i < 5 GeV, 
Piab < 12 GeV/c) in terms of interactions between known hadrons and their resonances. At higher energies, 
y/s > 5 GeV, the excitation of color strings and their subsequent fragmentation into hadrons are taking into 
account. 

The model is based on the covariant propagation of all hadrons considered on the (quasi-)particle level on 
classical trajectories in combination with stochastic binary scatterings, color string formation and resonance 
decay. It represents a Monte Carlo solution of a large set of coupled partial integro-differential equations 
for the time evolution of the various phase space densities of particle species i = N, A, A, etc. The main 
components of the model are cross sections of binary reactions, potentials and decay widths of resonances. 

The potential interaction is based on a non-relativistic density-dependent Skyrme-type equation of state 
with additional Yukawa- and Coulomb potentials. Momentum dependent potentials are not used. The 
Skyrme potential consists of a sum of two- and a three-body interaction terms. The two-body term, which has 
a linear density-dependence models the long range attractive component of the nucleon-nucleon interaction, 
whereas the three-body term with its quadratic density-dependence is responsible for the short range repulsive 
part of the interaction. The parameters of the components are connected with the nuclear equation of state. 
Only the hard equation of state has been implemented into the current UrQMD model. 

The impact parameter of a collision is sampled according to the quadratic measure (dW ~ bdb). At 
the given impact parameter, the centers of projectile and target are placed along the collision axis in such 
a manner that the distance between surfaces of the projectile and the target is equal to 3 fm. Momenta 
of nucleons are transformed in the system where the projectile and target have equal velocities directed in 
different directions of the axis. After that the time propagation starts. During the calculation each particle is 
checked at the beginning of each time step whether it will collide within that time step. A collision between 
two hadrons will occur if d < a tot /n, where d and a tot are the impact parameter of the hadrons and 
the total cross section of the two hadrons, respectively. After each binary collision or decay the outgoing 
particles are checked for further collisions within the respective time step. 

In the UrQMD model the total cross section a tot depends on the isospins of colliding particles, their 
flavor and the cm. energy. The total and elastic proton-proton and proton-neutron cross sections are well 
known |54j . Since their functional dependence on -^/s shows a complicated shape at low energies, UrQMD 
uses a table-lookup for those cross sections. The neutron-neutron cross section is treated as equal to the 
proton-proton cross section (isospin-symmetry) . In the high energy limit (y/s > 5 GeV) the CERN/HERA 
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parameterization for the proton-proton cross section is used |54j . 
Baryon resonances are produced in two different ways, namely 



i) hard production: N+N-> AN, AA, N*N, etc. 

ii) soft production: 7r~+p— > A , K~+p— > A*... 

The formation of s-channel resonances is fitted to measured data. Partial cross-sections are used to calculate 
the relative weights for the different channels. 

There are six channels of the excitation of non-strange resonances in the UrQMD model, namely NN — > 
7VA1232, NN*,NA*, A1232A1232, A1232N*, and A1232A*. The A1232 is explicitly listed, whereas higher 
excitations of the A resonance are denoted as A*. For each of these 6 channels specific assumptions were 
made with respect to the form of the matrix element, and the free parameters were adjusted to the available 
experimental data. 

Meson-baryon (MB) cross sections are dominated by the formation of s-channel resonances, i.e. the 
formation of a transient state of mass m = ^/shh, containing the total cm. energy of the two incoming 
hadrons. On the quark level such a process implies that a quark from the baryon annihilates an antiquark 
from the incoming meson. Below 2.2 GeV cm. energy intermediate resonance states get excited. At higher 
energies the quark-antiquark annihilation processes become less important. There, t-channel excitations of 
the hadrons dominate, where the exchange of mesons and Pomeron exchange determines the total cross- 
section of the MB interaction [55] . 

To describe the total meson-meson reaction cross sections, the additive quark model and the principle of 
detailed balance, which assumes the reversibility of the particle interactions are used. 

Resonance formation cross sections from the measured decay properties of the possible resonances up to 
cm. energies of 2.25 GeV for baryon resonance and 1.7 GeV in the case of MM and MB reactions have been 
calculated based on the principle. Above these energies collisions are modeled by the formation of s-channel 
string or, at higher energies (beginning at ^/s = 3 GeV), by one/two ^-channel strings. In the strangeness 
channel elastic collisions are possible for those meson-baryon combinations which are not able to form a 
resonance, while the creation of i-channel strings is always possible at sufficiently large energies. At high 
collision energies both cross sections become equal due to quark counting rules. 

A parameterization proposed by Koch and Dover [56 is used in UrQMD model for baryon-antibaryon 
annihilation cross section. It is assumed that the antiproton-neutron annihilation cross section is identically 
to the antiproton-proton annihilation cross section. 

The hadron-hadron interactions at high energies are simulated in 3 stages. According to the cross sections 
the type of interaction is defined: elastic, inelastic, antibaryon-baryon annihilation etc. In the case of inelastic 
collision with string excitation the kinematical characteristics of strings are determined. The strings between 
quark and diquark (antiquark) from the same hadron are produced. The strings have the continuous mass 
distribution f(M) oc 1/M with the masses M, limited by the total collision energy ^/s: M\ + M 2 < ^/s. The 
rest of the ^/s is equally distributed between the longitudinal momenta of two produced strings. 

The second stage of h-h interactions is connected with string fragmentation. The string break-up is 
treated iteratively: String — > hadron + smaller string. A quark-antiquark (or a diquark-antidiquark) pair 
is created and placed between leading constituent quark-antiquark (or diquark-quark) pair. Then a hadron 
is formed randomly on one of the end-points of the string. The quark content of the hadron determines its 
species and charge. In case of resonances the mass is determined according to the Breit-Wigner distribution. 
Finally, the energy-fraction of the string which is assigned to the newly created hadron is determined: After 
the hadron has been stochastically assigned a transverse momentum, the fraction of longitudinal momentum 
transferred from the string to the hadron is determined by the fragmentation function. The conservation 
laws are fulfilled. The diquark is permitted to convert into mesons via the breaking of the diquark link. 

This iterative fragmentation process is repeated until the remaining energy of the string gets too small 
for a further fragmentation. 

The fragmentation function f(x,mt) represents the probability distribution for a hadron with the trans- 
verse mass mt to acquire the longitudinal momentum fraction x from the fragmenting string. One of the 
most common fragmentation functions is the one used in the LUND model [57] . In UrQMD, different 
fragmentation functions are used for leading nucleons and newly produced particles, respectively: 



nuc 




for leading nucleons 



(3) 



/(^Oprod 



(1 — x) 2 for produced particles 



(4) 
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with A = 0.275 and B — 0.42. The fragmentation function /(x) pro d, used for newly produced particles, is 
the well-known Field-Feynman fragmentation function ~~~~J E~~~ ■ 

The fragmentation scheme determines the formation time of created hadrons. Though there are various 
possibilities (for details see Refs. [TT| "12]). 

After the fragmentation, decay of the resonances proceeds according to the branching ratios compiled by 
the Particle Data Group [!~~fl [37] . The resonance decay products have isotropical distributions in the rest 
frame of the resonance. If a resonance is among the outgoing particles, its mass must first be determined 
according to the Breit-Wigner mass-distribution. If the resonance decays into N > 2 particles, then the 
corresponding N— body phase space is used to calculate their N momenta stochastically. 

The Pauli principle is applied to hadronic collisions or decays by blocking the final state if the outgoing 
phase space is occupied. 

The UrQMD collision term contains 55 different baryon species (including nucleon, delta and hyperon 
resonances with masses up to 2.25 GeV) and 32 different meson species (including strange meson resonances), 
which are supplemented by their corresponding anti-particle and all isospin-projected states. The states can 
either be produced in string decays, s-channel collisions or resonance decays. For excitations with higher 
masses than 2 GeV a string picture is used. Full baryon/antibaryon symmetry includes: The number of the 
implemented baryons therefore defines the number of antibaryons in the model and the antibaryon-antibaryon 
interaction is defined via the baryon-baryon interaction cross sections. 

The elementary cross sections are fitted to available proton-proton or pion-proton data. Isospin sym- 
metry is used when possible in order to reduce the number of individual cross sections which have to be 
parameterized or tabulated. 

We have installed a bug fixed version of UrQMD 1.3 [33] at the Hep Web server. 



A role of heavy resonances in hadron-hadron interactions becomes more and more essential with a growth 
of a collision energy. The resonances have wide mass distributions. The distributions overlap, and the mass 
spectra in the binary reactions become practically continuous and smooth. This is taken into account in 
the Fritiof model ~~~~] . The model assumes that all hadron-hadron interactions at high energies are binary 
reactions: a + b — > a' + b' , where a' and b' are excited states of initial hadrons a and b. The excited states 
are characterized by masses which are selected by the following procedure: in the center of mass of colliding 
hadrons the energy-momentum conservation low has form: 



where E a and Ej, [E a i, Ey) ~ energies of initial (final) hadrons a and b (a',b r ), p az and Pf, z - longitudinal 
momentum components, p a >±_ and P b>± ~ transverse momentum components of hadrons a' and b' . 
Adding and subtracting the first two equations from (~~~| leads to: 
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where P + = E + p z , P 
At high energies: 
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The limits of P a 7 and P b t are defined as 



(9) 
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In case of hadron-nucleus interactions the kinematics governed by Eqs. |5]) - ([9]) is applied for the first 
collision of the projectile hadron with one of the target nucleons (a + Ni — > a' + N[). For the second collision 
(a' + N 2 — > a" + N 2 ) , analogous relations are used but ^ is changed by 

[Pa'^N 2 i [P^Pa'}- (10) 

As a result, the sequence of the collisions leads to a systematic increasing of the mass of the hadron a if 
transferred momenta are small. 

A similar approach is also applied to simulate nucleus-nucleus interactions. Here the reactions a' + b' — > 
a" + b" arc considered. The distributions on P~„ and P^, are the same as those for P~, and Pu , but the 
limits of P~„ and Py, are redefined as 

[P«,Py], [P£,P$]. (11) 

Thus, the model considers the interactions of the cascade particles between themselves from the point of 
the view of the cascade-evaporation model. Probabilities of the multiple interactions are calculated in the 
Glauber approach. 

The excited hadrons are treated as QCD-strings, and a corresponding algorithm is applied for a simulation 
of their fragmentation into observed hadrons. Due to the increasing of the masses, the multiplicities are 
increased too. Thus, the two factors - multiple interactions and the increasing of the masses, explain the 
increase of the produced particle multiplicity at a passage from hadron- nucleon interactions to the hadron- 
nucleus and nucleus- nucleus ones. 

In the version of the model code installed at the HepWeb server we take into account the elastic scattering 
of hadrons as well as the inelastic interactions (for details see Ref. |61j). 

The cascading of the produced hadrons in nuclei was not considered in the model. As a results, a 
multiplicity of slow particles associated with a destruction of residual nuclei was not described sufficiently 
well. In order to overcome the drawback it was proposed in Refs. [SHI HH] to enlarge the Frition model by 
the reggeon theory inspired model of the nuclear destruction (RTIM) [3T] . 

A simulation of the destruction consists from two stages. At the first stage a set of interacting nucleons 
(wounded nucleons) is determined with a help of the Glauber approximation. At the second stage non- 
interacting nuclear nucleons are considered. It is assumed that a non-interacting nucleon located at a 
relative impact parameter distance, r, from a wounded nucleon can be involved in the interactions with a 
probability 

W = C nd e- r2 / r ^, 

where C n d and r n d are parameters. The involved nucleon can involve another spectator nucleon, and so 
on. It is assumed also that all wounded and involved nucleons leave the nucleus. Good results have been 
obtained for light nuclei at C n d = 1 and r n d — 1.2 fm. We recommend to use C n d = 0.2 and r n d — 1 fm for 
heavy nuclei. 

To ascribe momenta of the escaped nucleon, we use the algorithm proposed in Ref. 29J . To explain its 
feature let us consider a reaction of a compound system, (1, 2), with a hadron h: (1, 2) + h — > 1 + 2 + h. 
Neglecting transverse momenta, a final state of the reaction will be fully determined by a value of merely 
one independent kincmatical variable. As the variable let us take 

x+ = {E 1 +p 1 )/(E 1 + E 2 + Pl +p 2 ). 

It is obvious that an analogous variable for the second particle, x 2 , will satisfy the condition - xf + x 2 = 1. 
All the other kinematical variables can be determined from the energy-momentum conservation law. 

In a case of a dissociation of two compound systems, A and B, containing A and B constituents, respec- 
tively, the z-th constituent of system A will be described by 

xf = (E Ai +Pi z )/WX and p i± , 
and the j-th constituent of system B - by 

yj = (E Bj - q ]z )/W B and q i± . 
Here, EAi(EB^) and Pi{qi) are energy and transverse momentum of i-th constituent of A(B). 

A B 
i=l i=l 
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Using the variables, let us write the conservation law as 
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\ - rnf ± 

h 4 
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(12) 



A 

E 

i=l 



Pi-L + E * 



where m?j_ = m 2 + pf t , = /if + and mi(Hi) is a mass of the i— th constituent of system A(B), 

The system (12 1 allows one to determine W A ,W B and kinematic properties of all the particles in the 
finite sets {x+,p~[}, {viiQi±}- 

W+ = (W ( ;W+ +a-P+ VA)/2Wq; (13) 
W B = (Wo W+-a + p + VA)/2W+ ; (14) 

W+ = (E A + E B ) + (P° Az +P Bz ); 

W - = (E A +E B )-(P° LZ +P Bz ); 

A 2 B 2 

« = E^. ^ = E^; 

i=l y « 



i=l 

\2 , J , o2 



= ^ - )A = -{W B y- - 

S W Vi W B 

According to experimental observations |62j . the average transverse momentum of a spectator fragment 
with mass number F obeys to the parabolic law: 

F(A-F) . 



< P'i >= ~ ^ ~ ' < P± >, y<Pi> = °- 05 GeV/c. 

To reproduce the result, the values of pi± for the knocked-out nucleons were sampled according to the 
distribution: 

dW oc exp{-fiJ < p 2 ± >)d 2 p l± , \j<p\ > = 0.05. (15) 

The sum of the transverse momenta (with "minus" sign) was ascribed to the residual nucleus. 
The choose of xf is carried out by 

dW cx exp[-(x+ - l/A) 2 /{d x /A) 2 ]dx+ , d x = 0.05. (16) 

The dispersion of the distribution was defined by fitting the average emission angle of 6-particles. x + of the 



+ 



residual nucleus was included as 1 — 

The knocked-out nucleons occurring in the zone of active action of nuclei was assumed to change again 
their characteristics. The new values of x~l and pi± were simulated using the distributions ( 15 ) and ( 16 ) at 
< p\ >= 0.3 (GeV/c) 2 and d x = 0.21. The results fr om |63| 1641 165] were used for a determination of the 
fitting parameters. 

The calculation of the nuclear residual excitation energy is carried out according to Ref. |66j . The 
implemented method is described in details in Ref. 29 . For a simulation of a relaxation of the excited 
nuclei the standard evaporation model is used [44]. 

A version of the model applied in Refs. [29, 67J underestimated negative charged particle multiplicities 
especially in the target fragmentation region. Thus, an attempt to take into account non-nucleon degree of 
freedom in nuclei has been considered. 

It is no doubt that the nucleons in nuclei can have virtual transitions - N — > N + ir and N + ir — > N. The 
virtual pairs can fall onto the mass-shell and become real ones due to an interaction with a projectile hadron. 
A calculation of cross sections and properties of the processes requires a solution of numerous questions of 
the nuclear theory, and can not be done in a full volume in the present time. Additionally it is needed 
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to consider a creation of the Nir pairs in the reggeon cascade. Thus, for a qualitative understanding of 
processes we have assumed that an involved or wounded nucleon can be the Ntt pair with a probability of 
~ 20%. More concrete, we assume that a proton or a neutron participated in the interaction can be A + - or 
A isobar. It is the most simple solution from a programming point of view. It throws away a complicated 
question on kinematical properties of the pairs. 



Appendix: HIJING and AMPT Models 

The quark-gluon degree of freedom manifests itself directly at super high energies in hadron-hadron inter- 
actions - jets of hadrons are produced which are interpreted as results of fragmentations of quarks and 
gluons suffered hard collisions with large transverse momenta. Inclusive properties of the jets are described 
by the quantum chromo-dynamics (QCD) with a high precision. The most popular Monte Carlo codes 
for a simalation of the processes are PYTHIA [5] and HERWIG [3]. Because the hadronisation process - 
the fragmentation of quarks and gluons into hadrons, can not be described in the frame of QCD, various 
phenomenological models are applied in the case. The total and elastic hadron-hadron cross sections are 
not described by QCD also. They are determined by the so-called soft interactions with low transverse 
momenta. Thus, for a complete interpretation of experimental data various models of a coupling of the 
hard and soft processes were proposed. One of the successful model is the HIJING code (Heavy Ion Jet 
Interaction Generator IT7] ). 

The model is a generalization of the Fritiof model. It considers the hard jet production simulated with 
a help of the Pythia code. It is assumed in the model that the inelastic nucleon-nucleon interaction cross 
section has a form: 

I _ e -(<T 80 /t+ow<i)Tjv(&) = fiY) 



Jn = d 2 b 



= d 2 b 



1 _ g-^.o/t T N (b) e -a hard T N (b) + 



[cThard T N (b)¥ ahard TN{ b) 



j=i J 

where Tjv is an overlap function of partons (quarks and gluons) from colliding hadrons at impact parameter 
b. It is given by the Fourier transform of the dipole form factor of the proton. 

T N (b,s) = 2^^, xo (0 = d(^KM), C = \b\Ms), 
(Jsoft(s) 96 

pi = 3.9, 7T &o = &soft(s)/2, (T so ft(s) = 57 (mb). 

K% - Mac-Donald function, ahard ~ integrated inclusive jet production cross section given by QCD. 

Eq. [17] allows one to determine a number of the hard scatterings accompanied all time by the excitation 
of hadron residues. If the number is equal to zero, the soft interaction is simulated as in the Fritiof model. 
For each hard interaction the kinetic variables of the scattered partons (gluons) are determined by calling 
PYTHIA 3] subroutines. Since jet production is dominated by gluon scatterings, it is assumed that quark 
scatterings only involve valence quarks, and hard processes are restricted only by the gluon-gluon scatterings. 
Simplification is also made for the color flow in the case of multiple jet production. Produced gluons are 
ordered in their rapidities and then connected with their parent valence quarks or diquarks to form string 
systems. Finally, fragmentation subroutine of JETSET is called for hadronization. 

Eq. [17] is entering the Glauber formulae for hadron-nucleus and nucleus-nucleus cross sections in the 
corresponding routines of the HIJING code. Besides a change of the gluon distribution in nuclei (the nuclear 
shadowing) and interactions of the scattered gluons with spectators nucleons (the final state interactions, or 
jet quenching) are taken into account. 

According to the Glauber approximation, it is assumed that a nucleus-nucleus collision can be decomposed 
into binary nucleon-nucleon collisions which generally involve the wounded nucleons. In a string picture, 
the wounded nucleons become strings excited along the beam direction. At high energy, the excited strings 
are assumed to interact again like the ordinary nucleon-nucleon collisions before they fragment. Unlike 
the Fritiof model, it is allowed that an excited string can be de-excited within the kinematic limits in the 
subsequent collisions. The binary approximation can also be applied to rare hard scatterings which involve 
only independent pairs of partons. The probability for a given parton to suffer multiple high scatterings 
is small and is not implemented in the current version of the program. A three-parameter Wood-Saxon 
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nuclear density is used for a calculation of the number of binary collisions at a given impact parameter. 
After all binary collisions, the scattered partons are connected with the corresponding valence quarks and 
diquarks to form string systems. The strings are then fragmented into particles. 

One of the most important nuclear effects in rclativistic heavy ion collisions is the nuclear modification 
of parton structure functions. It has been observed |68j that the effective number of quarks and antiquarks 
in a nucleus is depleted in the low region of x. This is called in the code as the nuclear shadowing, and is 
implemented. 

Another important nuclear effect on the jet production in heavy ion collisions is the final state integration. 
In high-energy heavy ion collisions, a dense hadronic or partonic matter must be produced in the central 
region. Because this matter can extend over a transverse dimension of at least Ra, jets with large Pr from 
hard scatterings have to traverse this hot environment. Thus it is very important to simulate the interaction 
of jets with the matter and the energy loss they suffer. It is estimated [BS] that the gluon bremsstrahlung 
induced by soft interaction dominate the energy loss mechanism. The induced radiation in HIGING is 
modeled via a simple collinear gluon splitting scheme with given energy loss dE/dz. It is assumed also that 
the interaction only occur with the locally comoving matter in the transverse direction. The interaction 
points are determined via a probability 

X s 

with given mean free path X s , where / is the distance the jet has traveled after its last interaction. The 
induced radiation is simulated by transferring a part of the jet energy AE(l) — IdE/z as a gluon kink to 
the other string which the jet interacts with. The procedure is continued until the jet is out of the whole 
excited system or when the jet energy is smaller than a cutoff below which a jet can not loss energy any 
more. The cutoff is the same as the cutoff Pq for jet production. To determine how many and which excited 
strings could interact with the jet, it is assumed that excited strings within a cylinder of radius r s along the 
jet direction could interact with the jet. X s and r s are two parameters in of the jet quenching effect. 

The code of the HIJING model operates quite fast. Thus it was widely used by RHIC and LHC experi- 
mental collaborations at designs of experimental setups. 

The AMPT model (A Multi-Phase Transport model) [JfJ, is built on the HIJING code. It considers 
additionally to the HIJING subjects the final state interactions between the partons. The interactions of 
the produced hadrons with nuclear matter is taken into account too. 

Appendix: Calculation of nucleus-nucleus interaction cross sections 
at high energy in the Glauber approach 

The amplitude of scattering a nucleus A on a nucleus B, when each of them transforms from the initial state 
\i) into the final states |/), is given in the Glauber's approach [7T],[72]- [75] by the expression 

. A B 

Fa,b(0) = 1 -J A d 2 be** i; (f A) f B \l-l[l[(l- 1 (b-s j + T k ))\i B ,iA), (18) 
n J j=ik=i 

where pa is the momentum of the projectile nucleus A, q is the transferred transverse momentum, b is 
the impact parameter, 7 is an amplitude of elastic NN scattering in the impact parameter representa- 
tion, {sj}, j = 1, 2, ...,A and {rfc}, k = 1,2,..., B are coordinates of nucleons within, respectively, A 
and B nuclei on the impact parameter plane. These coordinates are measured from the center of mass of 
each nucleus respectively, too. 



Starting from eq. 18 it is possible to find AB elastic scattering amplitude 

A B 



FAB(q)= l -7^ J d 2 be*^ {l-f[f[(l-l(b-s 3 +T k ))}- (19) 



j=lfc=l 

A B 

|2 



\ip A {f x,..., ? A )\ 2 \ii B {ii,...,i B )\ Y[d 3 rj Y[dH k , 

j=i k=i 



differential cross section 
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da el 

and total cross section 



-pel 
* AB 



(20) 



4-7T 

a tf = _ ImF eJ {Q y (21) 
Pa 

The cross section of the quasi-elastic scattering of the nucleus A when it is conserved, but other nucleus 
B undergoes all excitations including destruction too (^4 + B —> A + X) is given by the expression 



A B 



a(A + B^A + X)= fd 2 b{l-f[f[ (1 - 7(6 - g t + } ■ 

J i=i k=i 

A B 

{i-nn (1-1*0-^+^))}- 
3=1 k=i 

2 2 ~* -* i2 

■ • ■ ,»*a)I l^(ri,...r A )| |V>b(£i,---,£b)| • 

A A B 

i—l i—1 z—1 

Finally, the cross section of the production of new particles can be defined as 

A B 



(22) 



prod 
AB 



[ d 2 b{l-f[f[ (l-p(b-s l + T k ))}- (23) 

J i=lfe=l 

A B 

i^ A (fi,...,r A )i 2 \ii)B{ti,...,t B )\ 2 n d * ri 



p (b)= 1 (b)+ 1 *{b)- 1 {b) 1 *(b). 



Eq 23 can be re-written in the form where each of terms would be interpreted as a probability of some 
process 



^ i=l 3=1 

A B 



1 p(b-S t + f fc ) 5j +ffc) 

2 ' jS 1 - P(6 - «i + 7*) 1 - p(6 - si + r fe ) 

A B ~| 

■ n n ci--p(b-si+T m ))+...\ 

1=1 ra=l ) 
A B 

\i/j A (r u ...,r A )\ 2 yj B {tu...,t B )\ 2 Y[d 3 n Y[d 3 U. 



Here the first term in the first braces is interpreted as a probability that the only one inelastic collision 
between i-th nucleon from nucleus A and j-th nucleon from nucleus B takes place when all nucleons co- 
ordinates are fixed. The second term describes a probability of inelastic collision of the A;-th nucleon from 
nucleus B with i-th and j-th nucleons in A nucleus, etc. Nucleons involved in the collisions are called as 
"wounded" ones, but the others are called as "spectators" ones. 

To calculate all cross sections discussed above, it is necessary to have a function 7(6) and the square of 
the modulus of the ground state wave function \tp\ of A and B nuclei. 

The approximation 

TW = °NN ■ 4n B e ( 25 ) 
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is often used at E > 1 GeV/nucleon. Here (Jnn is the total cross section of NN interaction, a is the ratio 
of the real part to the imaginary part of the elastic scattering amplitude at zero momentum transfer, B is 
the slope parameter of the differential cross section of the elastic NN scattering. Expression 25 corresponds 
to the following parameterization of the elastic NN scattering amplitude in the momentum representation 

« = ^^(l-i«) e- Bf/2 - 

Sets of values ct^, a and B at various energies are presented in the compilations [35]. The Particle Data 
Group proposed the following parametrizations for the total and elastic NN cross-sections at high energies 

a(p)=A + Bp n + Cln 2 (p) + Dln(p), (26) 

where p is the laboratory momentum in GeV/c. The parameters - A, B, C, and D are given in Ref. |37| . 
Function \iPa\ is often given as 



\^Pa\ 2 = Y[p A (,n), (27) 



where p represents the one-particle density of a nucleus. In this case the aggregate of the nucleon coordinates 
does not meet self-evident demands 

A 

Vr1=0. (28) 



Taking into account this condition is named " an account of center of mass correlation" . 

If projectile is a nucleon, it is necessary to set a density pa as a delta-function pa(t) = S(r) in this case. 
The parameterization from Ref. |76j is used for deuteron 



3 

KM 2 = 5>e" ?2 /47 S (29) 
»=i 



71 = 225 (GeV/c)- 2 , ci = 0.178/(47r 7l ) 3/2 , 

72 = 45 (GeV/c)- 2 , c x = 0.287/(47r 72 ) 3 / 2 , 

73 = 25 (GeV/c)- 2 , c x - 0.535/(47r 73 ) 3/2 . 

Here r is a distance between a proton and a neutron within a deuteron. 
For 3 H, 3 He and 4 He nuclei \iPa\ 2 has been chosen as 

A 



iM^H^m—^^^ (30) 

i=l »=i [kKa) 

where Rzh = R 3 He = 1-81 /to, R^He = 1-37 /to. For all other nuclei (A > 6) a one-particle density has 
been defined as 

r-R A 

Pa(i~) = const/ (1 + e 5 ), (31) 

with Ra = 1.07 • A 1 / 3 /to, c = 0.545 /to. The center of mass correlation has been taken as in paper [77] . 

The calculation of the cross sections is performed in accordance with the method discussed in [IT] , where 
the statement 

A B 

»=i j=i 

is treated as a probability measure to find out different sets of nucleon coordinates in A or B nucleus. In 
this case, the cross sections are found as mean values over various sets of nucleon coordinates. Thus it can 
be written in the form 



tot 

° AB A 



— • £ / d 2 bET B (b,{rA},{tB}) (32) 

stat . 1 J 
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where 

A B 

E% = l-lH[(l- 7 (b-sl + f j )), 
i=ii=i 

and N sta t denotes a number of various sets of nucleon coordinates. Therefore, the accuracy of the calculation 



depends on the value of N sta t- (The expressions analogous to 32 can be written for all other cross sections.) 



The function E is called as a distribution over impact parameter for inelastic AB interactions: 



N sta 

\ s r^prod/ 



p (b) = — ■ > . Ksib, {rl}, {t B }), (33) 



where 



1=1 3 = 1 

In accordance with the expression |24| the algorithm for simulation of inelastic collisions consists of the 
following steps: 

1. Calculation and tabulation of the function P(b); 

2. Generation of the impact parameter b in accordance with P(b) distribution; 

3. Sampling of the coordinates of the nucleons within nuclei according to \^Ja\ 2 arL d I^bI 2 distribution 
functions; 

4. Sampling and storing pairs of nucleons interacting inelastically. For this purposes A ■ B random 
numbers £y (i = 1, . . . , A, j = 1, . . . , B), uniformly distributed in the [0,1.] interval, are chosen. If 
£ij < p(b — Si+fj) then it is considered that an inelastic collision takes place between i-th nucleon of 
A nucleus and j-th nucleon of B nucleus. 

The first step is performed only once with given nuclei A and B. The steps 2-4 are repeated as many 
times as it is necessary. 

The method described above has been realized as a set of FORTRAN routines. These routines and their 
assignments are presented in the table. 

NN name * Routines' function Subject 

"1 TOTAL * Calc. total cross sect. u% 

2 ELASTIC * Calc. elastic cross sect. a e i 

3 AB_TO_AX * Calc. cross sect, of reaction A + B -> A + X 

4 AB_TO_XB * Calc. cross sect, of reaction A + B -> X + B 

5 AB_TO_X * Calc. cross sect, of reaction A + B -> X 

* without process of produc- 
tion of particles 



* r< Q l„ „™oo c ar .+ ~t vrod 



6 PRODUCT * Calc. cross sect, of multi 



AB 



particle production 

7 GJXLSECT * Calc. all above cross sect. 

8 GL.STAR * Simulates inelastic events 

All of them can be called through the HepWeb server. 
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